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It is shown that high energy features of the cosmic ray spectrum, in particular the kink around 4 PeV and the corre- 
sponding change in spectral index, may be explained from interactions between highly energetic protons and relic Big 
Bang antineutrinos, if the latter have a rest mass of about 0.4 eV/c 2 . This explanation is supported by experimental 
data from extensive air-shower experiments. 
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PACS numbers: 13.85.Tp, 14.60. Pq, 98.70.Sa 



According to the Big Bang model of the evolving Uni- 
verse, large numbers of (electron) neutrinos and antineu- 
trinos have been around since the beginning of time. Dur- 
ing the very first second, these (anti-)neutrinos were in 
thermal equilibrium with the other particles that made 
up the primordial soup: photons, electrons, positrons 
and nucleons. Photon-photon interactions created e + e~ 
pairs, which annihilated into photon pairs. Interactions 
between (anti-)neutrinos and nucleons turned protons 
into neutrons and vice versa. 

This leptonic era came to an end when this thermal 
equilibrium ceased to exist. In the reaction 



P 



(1) 



the mass difference between the particles in the final and 
initial states amounts to 1.8 MeV. For this reaction to oc- 
cur, the energy available in the & e p center-of-mass system 
has to exceed this value. When the temperature dropped 
below 1 MeV, reaction (1) became energetically impossi- 
ble in most D e p encounters. In addition, because of the 
expansion of the Universe and the corresponding reduced 
density, the mean free path of the neutrinos became in- 
creasingly long, and soon the average neutrino lifetime 
exceeded the age of the Universe. The thermal equilib- 
rium between neutrinos and other particles thus came to 
an end, the neutrinos decoupled and the neutron/proton 
ratio froze out, determining the primordial He/H ratio. 

Since that moment, the wavelengths of the neutrinos 
have been expanding in proportion to the size of the Uni- 
verse. Their present spectrum is believed to be a Fermi- 
Dirac distribution characterized by a temperature of 1.95 
K Q . The present density of these Big Bang relics is es- 
timated at ~ 100 cm -3 for each neutrino flavor, nine 
orders of magnitude larger than the density of baryons 
in the Universe. 

In spite of this enormous density, relic neutrinos have 
until now escaped direct detection. They are acknowl- 
edged as a potentially very important component of dark 
matter, but this hinges on their mass, which is not known 
either. The recent observation of a zenith-angle depen- 
dence in the ratio of v e /vn induced interactions in the Su- 
perkamiokande detector is generally considered strong 



evidence for (at least one) non-zero neutrino restmass, 
and suggests a lower limit of the order of 0.03 eV. The 
upper limit on the v e massQ is 15 eV || and comes from 
endpoint measurements of tritium (3 decay and from the 
measured spread in arrival times of v^s from supernova 
SN1987a. 

The reason why relic neutrinos have yet to be observed 
is their extremely small kinetic energy, which makes it 
hard to find a process through which they might reveal 
themselves. In this Letter, such a process is proposed and 
its consequences are examined. 

Let us imagine a target made of relic P e 's and let us 
bombard this target with protons. Let us suppose that 
we can tune this imagined beam to arbitrarily high ener- 
gies. Then, at some point, the proton energy will exceed 
the value at which the pv & center-of-mass energy exceeds 
the combined mass of a neutron and a positron. Beyond 
that threshold, the inverse /3-decay reaction 



P- 



(2) 



is energetically possible. The threshold proton energy 
for this process depends on the mass of the D e target 
particles. If this mass is large compared to the kinetic 
energy of the target particles, this may be considered a 
stationary-target problem, and the center-of-mass energy 
can be written as 



E c .m — 



2E p m u 



2E p m u (3) 



When the experimental mass value of the proton (938.272 
MeV) is substituted in eq. [3] and the requirement is made 
that E ciri > m n + m e (940.077 MeV), this leads to 



E p m v > 1.7 • 10 ib (eV) 



(4) 



Thus, for a neutrino mass of 1 eV, this process may take 
place for proton energies in excess of 1.7 PeV, while the 
threshold energy is 17 PeV when the mass is only 0.1 eV. 



*In the following, we express masses in energy units, omit- 
ting the c -2 factor. 
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Obviously, the target needs to be huge to see any effect 
of the described process. The largest imaginable v e tar- 
get is the Universe. This target is continuously traversed 
by cosmic rays, 90% of which are free protons (the re- 
maining fraction consists predominantly of a particles). 
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FIG. 1. The all-particle spectrum of charged cosmic rays 

Figure [I] shows the spectrum of these cosmic rays, com- 
piled by the Particle Data Group [|| . This spectrum falls 
very steeply with energy. At energies below 1 PeV, the 
differential energy spectrum is well described by 



dN 
~a T E 



1.8 E 



-2.7 



nucleons 



cm 4 sr s 



GeV 



(5) 



In order to display characteristic features of this spec- 
trum which would otherwise be hard to discern, the dif- 
ferential energy spectrum has been multiplied by E 2 ' 7 
in this figure. The steepening that occurs between 1 
PeV and 10 PeV is known as the knee of the cosmic 
ray spectrum. Ever since its discovery, this feature has 
been the subject of intense interest and many explana- 
tions have been proposed. Especially models in which 
the cosmic rays are the result of particle acceleration in 
the shock waves produced in supernova explosions (the 
Fermi mechanism) have received much attention in the 
literature. Such models predict a maximum energy, pro- 
portional to the nuclear charge Z of the particles || . 

It is expected that the cosmic ray spectrum exhibits a 
high-energy cutoff at ~ 10 20 eV ||. Beyond that thresh- 
old, photopion production in collisions between cosmic 
protons and microwave background photons is energet- 
ically possible. Therefore, one does not not expect to 
observe any such protons, unless their source is located 
relatively nearby (i.e. , < 50 Mpc). 

The inverse /3-decay reaction (2) could, in principle, 
lead to a similar depletion of protons at energies above 
the threshold. If this process were indeed responsible for 
the change in the spectral index observed around 4 PeV, 
the neutrino mass should thus have a value of about 0.4 
eV (4). We therefore postulate the following hypothesis: 

The mass of the relic electron neutrinos is 0.4 eV. The 



change of the spectral index in the all-particle cosmic ray 
spectrum at an energy of ~ 4 • 10 15 eV is caused by the 



onset of the reaction p + v e — > n -\ 
energetically possible at this energy. 



e + , which becomes 



Let us examine this hypothesis. The most crucial ques- 
tion concerns the probability of the pv e reactions. In 
their historic experiment, Cowan and Reines measured 
the cross section for the inverse process (1), at energies 
just above the threshold, to be ~ 10~ 43 cm 2 J7J. As the 
proton energy increases above the threshold value, the 
cross section for the reaction p + D e — > n + e + rises pro- 
portional to the available energy in the center-of-mass 
system, i.e. the energy above the Q- value of 1.8 MeV [0: 

o- p p e ^ne+ = 0.34(£ cm -Q)10- 41 cm 2 (6) 

with the energies expressed in units of MeV. For proton 
energies between 10 16 eV and 10 18 eV, the cross section 
thus rises from 10~ 41 — 10~ 39 cm 2 . At a target density 
of 10 2 /cm 3 , this cross section translates for a 10 16 eV 
proton in a mean free path of 10 39 cm, or an average 
lifetime of 10 21 years, 11 orders of magnitude longer than 
the age of the Universe. Therefore, it seems very unlikely 
that the proposed process could affect the cosmic proton 
spectrum at a significant level. 

However, it is important to realize that, with a mass of 
0.4 eV, the relic Big Bang P e 's would be nonrelativistic 
(kT ~ 10~ 4 eV). Typical velocities would be of the order 
of 10 5 m/s in that case jjj, less than the escape veloc- 
ity from the surface of the Sun. Such neutrinos may be 
expected to have accumulated in gravitational potential 
wells. Weiler || recently estimated that the density of 
relic neutrinos in our own galaxy would increase by four 
orders of magnitude (compared to the universal density 
of 100/cm 3 ) if their mass were 1 eV. 

Locally, this effect could be much more spectacular. 
Extremely dense objects, such as neutron stars and black 
holes could accumulate very large numbers of relic neu- 
trinos and antineutrinos in their gravitational field. Let 
us consider, as an example, a typical neutron star, with 
a mass (M) of 3 • 10 30 kg and a radius of 10 km. Even at 
a distance (r) of one million kilometers from this object, 
the escape velocity is still considerably larger than the 
typical velocity of these relic neutrinos: 700 km/s. 

The concentration of relic neutrinos in such a local po- 
tential well is ultimately limited by the Pauli principle, 
which limits their phase-space density to 4g„/i~ 3 , where 
g v denotes the number of helicity states and h Planck's 
constant Q. Since the escape velocity scales with r -1 / 2 , 
the maximum neutrino density is proportional to r~ 3 / 2 , 
and reaches values of O(10 12 /cm 3 ) near the surface of 
this neutron star. If the source of the potential well has 
a larger mass, the achievable neutrino densities increase 
proportional to M 3 / 2 . In the "neutrino atmosphere" sur- 
rounding such dense objects, the lifetime of a 10 16 eV 
proton could thus be reduced to values comparable to, 
or even considerably smaller than the age of these ob- 
jects. 
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The example of a neutron star was chosen since objects 
of this type may well play a role in accelerating protons 
to the very high energies considered here. Neutron stars 
usually rotate very fast and exhibit very high magnetic 
fields (typically ~ 10 8 T). When the magnetic axis does 
not correspond to the rotation axis, the changing mag- 
netic fields in the space surrounding the neutron star may 
give rise to substantial electric fields, in which charged 
particles may be accelerated to high energies Q. The 
strong magnetic fields, in turn, will tend to confine the 
accelerated particles to the region surrounding the neu- 
tron star. When the charged particles are electrons, the 
acquired energy is emitted in the form of synchrotron 
radiation, a characteristic signature of many pulsars. 

However, for protons this type of energy loss only plays 
a significant role at energies that are 13 orders of mag- 
nitude ((7 P /7e) 4 ) larger than for comparable electrons. 
For example, a 10 16 eV proton describing a circular orbit 
with a radius of 1,000 km loses per orbit about 0.6% of 
its energy in the form of synchrotron radiation. To put 
this in perspective, that is about a factor of 5 less than 
the fractional energy loss experienced by 100 GeV elec- 
trons in CERN's LEP ring. It is, therefore, not inconceiv- 
able that a stationary situation might arise in which pro- 
tons and heavier ions, initially attracted by gravitational 
forces, and accelerated to high energies in the electric 
fields mentioned above, end up in a closed orbit around 
the neutron star in which energy gained and synchrotron 
losses balance each other. In that case, collisions with 
relic neutrinos accumulated in the gravitational field of 
the neutron star might become a realistic possibility. 

The sharpness of the knee suggests that, if our hypoth- 
esis is correct, most of the protons with energies above 
the threshold trapped in the gravitational field undergo 
this fate. One may then wonder why the cosmic ray spec- 
trum does not exhibit a "step" at the threshold energy. 
Two factors are responsible for this. 

First, there are the effects of the detector response func- 
tion. Upward fluctuations in the measured signal have a 
major impact on the measured spectrum of a very steeply 
falling distribution such as this one. Especially when this 
response function is not Gaussian, these effects may ex- 
tent over a considerable energy range. 

Second, the neutrons produced in (|^) eventually decay 
into protons which carry, on average, 1/6 of the energy 
of the protons that initiated the process, with a broad 
distribution about this average. Decay protons with an 
energy below 4 PeV contribute to a "bump" in the spec- 
trum just below this threshold. Such a bump is indeed 
experimentally observed, and is known as the kneecap 
(sec fig. [l]). As a matter of fact, this bump contributes 
significantly to the sharpness of the knee. 

As a result of these two effects, the cosmic ray spec- 
trum would be gradually depleted of protons at energies 
beyond the 4 PeV threshold. 

The cosmic ray spectrum is a rich source of possible 
tests of our hypothesis. For example, if the knee is indeed 



caused by the onset of the inverse /3-decay reaction, the 
observed kink should be a feature for protons and protons 
alone. But there is more. 

Other cosmic ray particles may interact with the relic 
background as well. The reactions in which a particles 
are dissociated in collisions with relic v e and v e 



a + v e — > 3p + n + e 
a + D e —> p + 3n + e + 



(7) 

(8) 



have Q-values of 27.5 MeV and 30.1 MeV, respectively. 
The threshold energies for these reactions are larger than 
the threshold energy for reaction (g) by factors of 61 and 
66, respectively. 

The cosmic ray spectrum in the range around the knee 
was measured in detail by the Akeno extensive air-shower 
experiment PH| . From a fit to their experimental data, 
the kink was located at \ogE p (eV) = 15.67. If our ex- 
planation of the origin of this kink is correct, a kink in 
the spectrum of cosmic a particles should thus be ex- 
pected at an energy that is larger by a factor of 61-66: 
logi^Q w 17.5 (E a « 3 • 10 17 eV). By contrast, supernova 
Shockwave models predict kinks in the proton and he- 
lium spectra separated by only a factor of two (Zn e /Z-a) 
in energy. Accurate measurements of the p/a ratio in the 
energy range 10 14 — 10 18 eV should thus make it possible 
to obtain conclusive evidence either way. 

It is remarkable that several extensive air-shower ex- 
periments have reported a kink around log-E = 17.5. 
The Fly's Eye detector, which obtained the largest event 
statistics, observed a change in the spectral index from 
3.01 ± 0.06 for E < 10 176 eV to 3.27 ± 0.02 for energies 
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FIG. 2. Experimental results obtained with the Fly's Eye 
detector \121. See text for details. 
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10 i7.6 < E < 10 i8.5 eV ||_ The Haverah Park detector 
also reported a kink at log-E = 17.6, with the spectral 
index changing from 3.01 ± 0.02 to 3.24 ± 0.07 |j. 

The cosmic ray spectrum obtained with the Fly's Eye 
detector is shown in fig. |^6. This detector is also ca- 
pable of measuring the altitude (and thus the amount of 
air traversed by the shower particles, A max ) at which the 
shower maximum is located (fig. |^a). This distribution 
exhibits an unmistakable kink at approximately the same 
energy as the kink in the cosmic ray spectrum, indicating 
an abrupt change in chemical composition. The dashed 
line in fig. ||a represents the results of a simulation in 
which the authors modeled the showers initiated by pro- 
tons [[l2| . To estimate what the equivalence of the dashed 
line will look like for showers induced by a particles, two 
effects have to be considered. 

First, the interaction length (A; nt ) of a's is shorter than 
for protons. For protons in air, Aj n t amounts to 90 g/cm 2 . 
Since the total cross section for interactions induced by 
highly energetic ions with atomic number A scales with 
A 2 / 3 , A int scales with A -2 / 3 , which gives A int = 55 g/cm 2 
for a's in air, a difference of 35 g/cm 2 . 

Second, the energy carried by the incoming particle is 
distributed to a larger number of secondaries. This has 
the effect that, in the subsequent shower development, 
the shower maximum is reached earlier than for proton 
showers of the same energy. The effect of this can, in 
first approximation, be estimated by considering the ion- 
induced shower as a collection of A showers, each with an 
energy E^/A. From the slope of the dashed line in fig. 
^|a, one may estimate a induced showers in air to reach 
their maximum earlier than proton induced ones, by 30 
g/cm 2 , counting from the first nuclear interaction. 

These (somewhat oversimplified) considerations thus 
lead to A max values for a induced reactions in the atmo- 
sphere that are 65 g/cm 2 smaller than those for protons 
of the same energies. This result is represented by the 
dotted curve in fig. ^p,. 

If the hypothesis described in this Letter were correct, 
one would expect the X max distribution to exhibit the 
following features. At energies below 10 15 eV, the X max 
value follows a line parallel to the ones for protons and 
a's, somewhere in between these two, depending on the 
relative H/He abundances. At 4 • 10 15 eV, the threshold 
for reaction (2) is crossed. Beyond that point, protons 
gradually disappear from the spectrum. The X max curve 
thus shows a kink at this point, beyond which the slope 
is less steep than before. For energies around 10 17 eV, 
almost all protons have disappeared and the spectrum is 
dominated by a particles. At ~ 3 • 10 17 eV, a new thresh- 
old is crossed, this time for reactions (7) and (8). Beyond 
that point, a's gradually disappear from the spectrum. 
In that process, protons carrying a significant fraction of 
the energy are produced. This leads to a gradual enrich- 
ment in protons (i.e. , to a slope in the X max distribution 
that is steeper than that for protons or a's) . Beyond 
10 18 5 eV, protons of extragalactic origin add to the mix. 

The Fly's Eye's experimental data (fig. ||a) are in 



detailed agreement with the described scenario. Fur- 
ther measurements of A max , at energies in the range 
10 15 — 10 17 eV, would constitute another important test. 



If the hypothesis decribed in this Letter turned out to 
be correct, a new era in cosmic ray physics may be en- 
visaged. Studies of the multi-PeV proton component of 
the spectrum would give a direct window on a very early 
epoch in the history of the Universe, i.e. the leptonic 
era. Proof of the existence of relic neutrinos would be 
a very strong additional argument in favor of the Big 
Bang model. Further measurements of their properties 
could provide more detailed insight in the processes that 
shaped the Universe in its first second. A neutrino mass 
of 0.4 eV would crucially increase the understanding of 
dark-matter issues. A 0.4 eV neutrino mass would also 
provide new perspectives for neutrino astronomy. For ex- 
ample, if a supernova occurred at a distance of 1 Mpc, the 
arrival times of simultaneously produced neutrinos of 1 
MeV and 10 MeV would differ by about 7 seconds. Neu- 
trino astronomy might one day even provide a precision 
tool for distance measurements on this basis. 

The high-energy cosmic proton and a spectra may thus 
hold the key to an important breakthrough in our under- 
standing of the physical world. 
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